Odd-frequency triplet pairing in mixed-parity superconductors 
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We show that mixed-parity superconductors may exhibit equal-spin pair correlations that are 
odd-in-time and can be tuned by means of an applied field. The direction and the amplitude 
of the pair correlator in the spin space turn out to be strongly dependent on the symmetry of 
the order parameter, and thus provide a tool to identify different types of singlet-triplet mixed 
configurations. We find that odd-in-time spin-polarized pair correlations can be generated without 
magnetic inhomogeneities in superconducting/ferromagnetic hybrids when parity mixing is induced 
at the interface. 



Introduction Symmetry breaking is a central concept in 
physics, superconductivity being an exemplary case. In 
a conventional superconductor gauge symmetry is bro- 
ken and Cooper pairs condense in the most symmetric 
configuration. The discovery of superfluid '^He started 
the era of unconventional superconductivity where other 
symmetries get broken and alternative pairing mecha- 
nisms emerge with respect to the phonon mediated one. 
This happens in high-temperature cuprates, f-electrons 
superconductors, cobaltates, ruthenates, pnictides, or- 
ganics and many other superconductors [ll, [2| ■ 

The issue of the gap symmetry and the pairing mecha- 
nism has become more intriguing since the recent discov- 
ery of superconductivity in a number of materials whose 
lattices lack inversion symmetry with important im- 
plications for the symmetry of the superconducting state. 
In this framework, difficulties may emerge in detecting 
parity mixing of Cooper pairs partly due to the fact that 
conventional experimental approaches utilized for the de- 
termination of parity of Cooper pairs, such as Knight 
shift measurements, do not provide useful information 
when spin-orbit coupling is larger than the pairing en- 
ergy scale. 

In some of the unconventional superconductors time 
reversal symmetry can also be broken both for spin sin- 
glets and spin triplets as well as in the cases of supercon- 
ductivity coexisting with ferromagnetic long-range order 
Concerning this symmetry, in the majority of the su- 
perconductors the order parameter is assumed to be even 
in the frequency domain, so that it may be even or odd 
in space depending on whether the Cooper pairs form 
spin singlets or triplets, respectively. However, if we al- 
low for an odd-in-time dependence, more exotic types of 
pairings are possible, like that originally predicted in the 
context of liquid ^He For these states the triplet pair 
correlation function is symmetric under the exchange of 
spatial and spin coordinates but antisymmetric under the 
exchange of time coordinates. 

A significant advance in the understanding of odd-in- 
time pairing has been achieved since many experimental 
observations have been accounted for as a manifes- 



tation of odd-in-time equal-spin (OTES) pair correlations 
generated at the interface between a spin singlet super- 
conductor and a ferromagnet whose magnetization has 
an inhomogeneous profile both in amplitude and phase 
Him. The interplay of magnetic non-collinearity and 



odd-in time pairing emerges also in magnetically active 
interfaces that are able to change the spin direction of 
incident electrons across the junction [l^, [l3|- Other- 
wise, zero-spin projection odd-in-time pair correlations 
take place in normal/superconductor junctions between 
triplet superconductors and diffusive normal metals 14 1 
or ferromagnets [l^ as well as in conventional ballistic 
hybrids without spin-triplet ordering [lit . 

The main interest in the search for odd-in-time pair 
correlations is motivated by the possibility of engineering 
spin-polarized supercurrents travelling through strong 
ferromagnets on distances that do not depend on the 



strength of the exchange field [17|. The issue of func- 
tionalizing spin-polarized pair currents in superconduct- 
ing/ ferromagnetic hybrids is indeed a challenging task for 
the high potential towards innovative applications in the 
field of the spintronics [3| where it is the spin degree of 
freedom rather than the charge one to be exploited. 

The aim of this letter is to demonstrate how odd- 
in-time equal-spin pair correlations are intimately con- 
nected to the symmetry breaking in unconventional su- 
perconductors. We show that OTES pair correlations 
manifest in mixed-parity superconductors and their ex- 
istence is strongly dependent on the direction of an ap- 
plied field in such a way to single out different types 
of parity mixing configurations. We also find a signif- 
icant connection with the generation of spin-polarized 
pairs in hybrids. A different mechanism, with respect 
to those mentioned above, is presented for the case of a 
superconductor/Rashba metal/ferromagnet hybrid with- 
out requiring any inhomogeneity in the magnetic profile. 
Indeed, if parity mixing is yielded at the interface, i.e. via 
antisymmetric spin-orbit interaction, then OTES corre- 
lations are not trivial in all directions of the spin space 
assuming that the magnetization and the vector defining 
the spin-orbit interaction are not orthogonal. 
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Model and results We start by considering a generic 
superconducting system having an equal time order pa- 
rameter, i.e. even in time, for both the spin singlet dok 
and the triplet channel dk = {dik, c?2k, c?3k} expressed 
through the following matrix Ak: 

'-rfik + «d2k dsk + rfok 

where c?ok = c^o-k and dk = — d^k 0, HOl- For the 
classification purposes it is useful to introduce the vec- 
tor = rk -I- Pk, with Tk = dok -I- dok* dk and 
P]^ = idk X dj^, that, for a non-vanishing amplitude, 
yields a non-unitary superconducting state with two non 
equivalent branches in the excitation spectrum. This is a 
direct consequence of the product Ak Aj^. being not pro- 
portional to the identity matrix. The q-vector is decom- 
posed in two components, one lying in the plane defined 
by the vectors [dk, dk*] and the other perpendicular to 
it. The in-plane term identifies a non- unitary state with 
pure mixed singlet-triplet character, while the perpendic- 
ular component is proportional to the magnetic moment 
of the Cooper pairs. 

To describe such type of superconducting states in sys- 
tems that can also lack the inversion symmetry one may 



use the following generalized BCS model: 
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FIG. 1. (Color online) Schematic representation of the non- 
zero components of the odd-in-time equal-spin pair correla- 
tion function for different types of mixed parity configura- 
tions, a) and b) stand for the symmetric (oval) and asym- 
metric (trapezoid) spectral components for the two spin pro- 
jections of the odd-in-time equal-spin correlations along a 
given direction, c) and d) indicate direction and amplitude of 
the odd-in-time equal-spin correlations for a centrosymmetric 
mixed parity unitary and non-unitary superconductor with 
respect to an applied magnetic field, respectively, e) and f) 
depict direction and amplitude of the odd-in-time equal-spin 
correlations for a non-centrosymmetric mixed parity without 
and with and applied field, respectively. The light rectangle 
indicates the plane perpendicular to h in c) and d), while it 
is the plane where and dk lie in e) and f). See the main 
text for the definition of rk,gkirfok, and dk. 



= ^ ('fkCro -I- .gk • (Ts.s' + h • t7s,s')CksCks' + 
k,s,s' 

^(Ak,ss'CksC-ks' A^^s^clkscL') (1) 

k,ss' 

where uq is the identity matrix, a the Pauli matrices, 
Ck = £k ~ M is the kinetic energy measured with respect 
to the chemical potential /i, Ak,ss' £^re the equal time 
components of the order parameter and h is the mag- 
netic field, respectively, gk is the antisymmetric spin- 
orbit term due to Rashba interaction associated with the 
breaking of the parity symmetry and such that gk = 
— 5-k due to time-reversal invariance. For the quantita- 
tive analysis we assume a specific tight-binding spectrum 
in two dimensions of the form Ek = —4 [cos(A:2;)-f cos(A:j,)], 
in units of the hopping amplitude, though the results do 
not depend on the details of the free electron dispersion 
and on the dimensionality. 

To search for OTES we explicitly determine for the 
model Hamiltonian in Eq. [T]the odd- in-time local corre- 
lator for equal-spin pairing that reads: 



(2) 



where V is the volume system, Tt is the time ordering 
operator, (.?.) the thermal average and Ui = t,i along 
the direction i = x,y,z in the spin space, respectively. 

The results obtained are reported in the Table |T] where 
the conditions and the directions of non zero OTES are 
indicated. It is worth pointing out that even-in-time pair 
correlations are also present for all the cases associated 
with equal time non-trivial order parameters and they 
can be as well induced by an applied field. 

Let us firstly consider the cases of centrosymmet- 
ric superconductors. In order to understand how odd- 
frequency pair correlations occur in parity mixing it is 
useful to start with the singlet and triplet states sep- 
arately. For the pure spin singlet superconductor the 
OTES correlations occur only in the presence of an ap- 
plied field h and with spin projections in the plane per- 
pendicular to it. This can be understood considering 
that the effect of the field along a given direction, i.e. z, 
is to induce odd-in-time correlations in the Sz = triplet 
channel of the Cooper pairs 1^, 12|. Then, due to the 
SU{2) spin symmetry, they lead to non-zero amplitude 
for the equal-spin correlators at a generic angle 9 in the 
plane perpendicular to h. Due to the symmetry of the 
ground and the excited states the correlation function 
is symmetric for Z2 inversion along a given direction in 
the spin space (see Fig. On the other hand, if we 

consider the case of a pure spin triplet superconductor, 
though the spectral amplitude of the time dependent cor- 
relator is non-trivial for a given k-point, due to its odd in 
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TABLE I. Cases with non-trivial odd-in-time equal-spin triplet pair correlations. Y and N indicates non-zero and null amplitude 
of the corresponding variable, respectively, do and dk indicate the spin singlet and triplet part of the order parameter, rk and 
gk are the vectors associated with the non-unitary parity mixing component of the Qk-vector and the Rashba interaction in 
non-centrosymmetric system, respectively, h is the applied field. F" is the odd-in-time correlator for pairs with equal-spin a 
along a given direction. Zg indicates the directions coplanar to the gk-vector. 
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space parity the average over the Brillouin zone leads to 
a cancellation of the odd-in-time pair correlations. This 
result does not depend on the unitary character of the 
triplet superconductor. 

Taking into account the previous analysis, one can di- 
rectly construct the results for centrosymmetric mixed 
parity superconductors. There are two relevant cases to 
consider with respect to the unitary nature of the ground 
state, corresponding to a non-zero amplitude of the vec- 
tors r and p. If we have a unitary superconductor, OTES 
correlations occur in the presence of an applied field and, 
as for the case of the pure singlet, they are Z2 symmetric 
and in the plane perpendicular to h (see FiglT];). More 
interesting is the case of a non-unitary mixed parity state. 
For such configurations the main role is played by a non- 
zero amplitude of the vector r. This vector is non-trivial 
if time-reversal is broken for either the triplet or the sin- 
glet component. Indeed, for the spin triplet it can lack 
time reversal invariance either in the orbital or in the spin 
channel. The latter would imply also a non-zero ampli- 
tude for the p component. Alternatively, non- unitary r 
states can also be related to a non-trivial phase relation 
between the singlet and the triplet order parameter that 
cannot be eliminated via a global gauge transformation. 
For these conditions, as schematically depicted in Fig. 
[TJi, if the applied field is not perpendicular to r the OTES 
occur in all directions (parallel and perpendicular to h) 
and they also break the Z2 symmetry along a given spin 
projection axis (Fig. [TJj). Depending on the material sys- 
tem and on the k-symmetry of the order parameter, if it 
is possible to choose a magnetic field such as h _L ta, for 
all k then there will be a field direction for which OTES 
will lie in a plane otherwise they have a finite component 
for any space direction. It is worth pointing out that the 
angular dependence of OTES correlations versus h can 
be used to establish whether mixed parity states exhibit 
time reversal symmetry breaking and get insight about 
the k-structure of the gap amplitude. A possible ap- 



plication is provided by the system Pr(Osi_a:Rua;)4Sbi2 
where a mixed parity phase that violates time reversal 
symmetry 21j has been proposed to interpolate between 
the singlet and the triplet configurations achieved at the 
points X = and x = 1, respectively. As a general re- 
mark, the observation of OTES can be an alternative 
method to muon spin rotation experiments for the detec- 
tion of time reversal symmetry breaking states in mixed 
parity configuration. 

Let us consider the case of noncentrosymmetric super- 
conductors, thus mixed-parity states in the presence of 
an antisymmetric Rashba type interaction. In this re- 
spect, as pointed out in Ref. 24], the most convenient 



configuration for the triplet dk vector is such that dk| |gk, 
though a possible misalignment between the two vectors 
is not completely prohibited. The analysis has been per- 
formed in all the possible configurations ^231] while the 
present results, due to the limited space, refer only to 
the physical relevant ones. In Fig[T]3, we have sketched 
the zero-field distribution of OTES correlations in the 
absence of an external source of time reversal symmetry 
breaking. In this case, OTES are coplanar to g^. and 
exhibit a Z2 symmetric behaviour at any given spin di- 
rection. The application of a magnetic field modifies the 
distribution of OTES pair correlations in the spin space 
in a way that depends on the amplitude of the prod- 
uct h • g. Only when the field has a parallel component 
with respect to g then OTES pair correlations occur both 
along the direction of h and perpendicular to it. A simi- 
lar behaviour can also happen if the noncentrosymmetric 
system is non- unitary (see Fig. [If). In this case, the non- 
orthogonality of the field with respect to r also leads to 
OTES pair correlations in all the spin directions. 

A relevant aspect for getting OTES in all spin direc- 
tions is related to the structure of the excitations spec- 
tra. Indeed, the spectra at each k-point in the Brillouin 
zone can have four eigenvalues that fulfill or lack the 
mirror symmetry when turning E into —E. The most 
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general conditions for which ffk, i-e. the k-component 
of the H matrix, has two couples of opposite eigenval- 
ues are given by a vanishing value for the linear and 
the cubic coefficients of the characteristic polynomial 
Pfc = Det[H]^~Xl4] = aok+aikA + a2kA^-|-a3kA^-|-a4kA'', 
with ask = —TrH]f^ = and aik — — 4(h-rk + 2^fch-gk + 
Pk ' Sk)- If "^sk and oik are both zero, one can exactly 
demonstrate, by means of symmetry arguments only, 
that the longitudinal component of ¥"^^{1) {F^^^it)) 
with respect to the field (g- vector) is identically zero [23| . 

It is thus the removal of mirror symmetry in the en- 
ergy spectrum together with the parity of the OTES k- 
component to yield missing cancellation for the longitu- 
dinal part with respect to the field when averaging in 
momentum space. 

Conclusions We have shown that OTES pair correla- 
tions can occur in superconductors with parity mixing. 
For centrosymmetric systems this requires the presence 
of an applied magnetic field and can lead to an asymmet- 
ric distribution for opposite spin projections depending 
on the field orientation with respect to the mixed parity 
non- unitary component. In materials that lack inversion 
symmetry, OTES can occur in zero field with components 
coplanar to the gk-vector or along all spin directions 
when a magnetic field not orthogonal to the gk-vector 
is applied. Due to the intrinsic differences of the OTES 
correlations depending on the magnetic field direction, it 
is possible to distinguish between various types of mixed 
parity superconducting states. We find that the removal 
of quasi-particle degeneracy under parity and time sym- 
metry is a key aspect for avoiding cancellation in the 
momentum space and get OTES pair correlations in all 
spin directions. 

We point out that a possible way to detect OTES is 
by setting up an experiment where spin-polarized long 
range Josephson current is observed across a diffusive 
half-metallic ferromagnet at different uniform orienta- 
tions of the magnetization. For the cases discussed in 
the paper there is no need to have an inhomogeneous 
profile of the magnetic pattern or spin active interfaces. 
The direction of the magnetization is used as an effective 
field to generate OTES along a given direction and at 
the same time as a filter to select the OTES pair corre- 
lations along a given spin direction. In particular for the 
case where OTES are longitudinal to the applied field it 
is possible to have spin-polarized pairs crossing the fer- 
romagnet without the need of inhomogeneous profile in 
the magnetization. 

This result can lead to a significant advance in the 
context of hybrid systems where the occurrence of mixed 
parity and subdominant components close to the inter- 
face can be the sources for time and parity symmetry 
breaking in the energy spectra and in turn for equal-spin 
odd-frequency pair correlations even without magnetic 



inhomogeneities. This observation can be used to en- 
gineer solutions where long range pairing currents can 
be induced by including sources of parity mixing at the 
interface between singlet superconductors and ferromag- 
nets, especially in the form of Rashba type junctions, and 
by properly selecting the magnetization direction. Since 
this parity mixing mechanism works for a uniform mag- 
netic profile it can also account for the persistence of the 
long range supercurrents in hybrids even when the ap- 
plied magnetic field is such that no magnetic domains 
are present in the ferromagnet. 
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